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Abstract 

We present the first measurement of decay amphtudes in B (pK* and measurements of branch- 
ing fractions in i? — > (pK^*^ decays based on 78.1 fb~^ of data recorded at the T(45') resonance with 
the Belle detector at the KEKB e+e" storage ring. The decay amplitudes for the different <^if*° 
helicity states are measured from the angular distributions of final state particles in the transver- 
sity basis. The longitudinal and transverse complex amplitudes are |^oP = 0.43 ib 0.09 it 0.04, 
|A_l|2 = 0.41 ± 0.10 ± 0.04, arg(^ii) = -2.57 ± 0.39 ± 0.09, and arg{A±) = 0.48 it 0.32 it 0.06. The 
direct CP- violating asymmetries are found to be consistent with zero. 

PACS numbers: 13.25.Hw,14.40.Nd 
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B meson decays involving b sss transitions, such as B ^ (pK and (t)K* ^ are forbidden 
to first order in the Standard Model (SM), but proceed by second order loop diagrams 
(penguin and box diagrams), which lead to the fiavor changing neutral current transition 
b s. These processes provide information on the Cabibbo-Kobayashi-Maskawa matrix 
element Vts 01 and are sensitive to physics beyond the SM such as R-parity violating SUSY 
contributions to 6 ^ sss They can also be used to perform independent measurements 
of the CP- violating parameter sin20i Jslj. The branching fractions of B ^ (pK have been 
predicted by QCD-factorization |4j and PQCD The decay B — > (pK* is a mixture of 
CP-even and CP-odd states; polarization measurements allow us to project out the different 
CP states statistically. 

In this letter, we report the first measurement of the helicity state amplitudes in B^ — >■ 
(pK*^ decay by a full three-dimensional angular analysis. We also report measurements of 
branching fractions and direct CP asymmetries in B~^ —>■ (pK^, B^ <pK^, 
and B^ (pK*^ decays (charge conjugate modes are included everywhere unless otherwise 
specified). 

This analysis is based on a data set with an integrated luminosity of 78.1 fb~^ taken at 
the T(4S') resonance recorded by the Belle detector p at the KEKB e^e~ collider 0. This 
luminosity corresponds to (85.5 ± 0.5) x 10^ produced BB pairs. The beam energies are 8 
GeV for e~ and 3.5 GeV for e+. 

The Belle detector is a general purpose magnetic spectrometer equipped with a 1.5 T su- 
perconducting solenoid magnet. Charged tracks are reconstructed in a Central Drift Cham- 
ber (CDC) and a Silicon Vertex Detector (SVD). Photons and electrons are identified using 
a CsI(Tl) Electromagnetic Calorimeter (ECL) located inside the magnet coil. Charged par- 
ticles are identified using measured dE/dx in the CDC as well as information from Aerogel 
Cherenkov Counters (ACC) and Time of Flight Counters (TOP). A kaon likelihood ratio, 
P{K/n) = Ck/{Ck + Ct,), has values between (likely to be a pion) and 1 (likely to be a 
kaon), where C.k{tt) is derived from dE/dx, ACC and TOP measurements. 

Candidate (p —>■ K^K~ decays are found by selecting pairs of oppositely charged tracks 
that are not pion-like {P{K/tt) > 0.1). The vertex of the candidate charged tracks is 
required to be consistent with the interaction point (IP) to suppress poorly measured tracks. 
In addition, candidates are required to have a K^K~ invariant mass that is less than 10 
MeV/c^ from the nominal (p meson mass. 

Pairs of oppositely charged tracks are used to reconstruct — >■ tt^tt^ decays. The tt+tt" 
vertex is required to be displaced from the IP by a minimum transverse distance of 0.22 cm for 
high momentum (> 1.5 GeV/c) candidates and 0.08 cm for those with momentum less than 
1.5 GeV/c. The direction of the pion pair momentum must agree with the direction defined 
by the IP and the vertex displacement within 0.03 rad for high-momentum candidates, and 
within 0.1 rad for the remaining candidates. 

Charged tracks with P{K/tt) > 0.4 (< 0.9) are considered to be kaons (pions). For 
7r° — > 77, a minimum photon energy of 50 MeV is required and the 77 invariant mass must 
be less than 16 MeV/c^ from the nominal 7r° mass. K* candidates are reconstructed in three 
decay modes: K+n', K*+ K+tt^ and ir^7r+. The invariant mass of the 

K* candidate is required to be less than 70 MeV/c^ from the nominal K* mass. 

A B meson is reconstructed from a (p meson candidate and a, K 01 K* candidate and 
identified by the energy difference AE = E]^ — Pb^^, and the beam constrained mass 

Mbc = (PbJ^m)^ ~ {p'j^^y. Pbeam the beam energy in the center-of-mass system (cms) 
of the T(45') resonance, and E'^^ and p^*^ are the cms energy and momentum of the 
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reconstructed B candidate. The S-meson signal window is defined as 5.27 GeV/c^ < Mbc < 
5.29 GeV/c^ and |AE| < 64 (60) MeV for B ^ (pK {B ^ (pK*). The signal window is 
enlarged to -100 MeV < AE < 80 MeV for 5+ ^ ^K*+{K*+ K+n^) because of the 
impact of shower leakage on AE resolution. An additional requirement cos 6k* < 0.8 is 
applied to reduce low momentum 7r° background, where 6^* is the angle between the K* 
direction and its daughter kaon defined in the K* rest frame. In the signal window about 
1% of the events have multiple candidates. We choose the best candidate according to the 
value of the B vertex ■ 

The dominant background is continuum e+e~ — >• qq production (g = u,d,c,s). Several 
variables are used to exploit the differences between the event shapes for continuum qq 
production (jet-like) and for B decay (spherical) in the cms frame of the T(4S') [3]. These 
variables are combined into a single likelihood ratio Rs = Cs/{Cs + Cqq)^ where Cg {^qq) 
denotes the signal (continuum) likelihood. An additional variable cos 6h, which is the angle 
between the momentum and the daughter kaon momentum in the rest frame, is included 
for the (j)K~^ and (pKg channels. 

Backgrounds from other B decay modes such as B ^ KKK^*\ B — ^ /o(980)ii'*^*^(/o 
K~^K~), B (pKn, B — >■ KKKn, and feed-across between (pK* and (pK decay channels 
are studied. The contributions from B KKK^*^ and B /o(980)K(*)(/o K+K') are 
estimated from the K^K~ invariant mass distribution. The K^K~ mass distribution for 
B — > KKK^*^ is determined by Monte Carlo (MC) simulation assuming three-body phase 
space decay. The shape for /o(980) is obtained from MC, where a Breit-Wigner with a 40 
MeV/c^ intrinsic width is assumed. These contributions are estimated separately by fits to 
the events outside of the (p mass region. The contribution from B KKK^*^ is estimated to 
be 5 — 9% of the signal yield and is subtracted from the raw sig nal yield. The B f^K^*^ 
contribution is estimated to be 2 — 12%. The large uncertainty in the intrinsic width of the 
/o(980) is included in the systematic error. The background from B —>■ cpKir decay, as well 
as higher K* resonance decay, is studied by performing fits to the Kt^ invariant mass. The 
estimated background (1 — 3%) is considered as a systematic error. Contamination from 
four-body B — > KKKir decays is checked by performing fits to the non-resonant region of 
K^K~ and Kit mass. It is found to be very small and is neglected. The feed-across from 
(pK* in (pK decay is removed by excluding events with AE < —120 MeV from the fit. A 
veto is applied in i? — > (pK* channels to remove the feed-across background from cpK. 

The signal yields {Ng) are extracted by extended unbinned maximum-likelihood fits per- 
formed in AE and M^c simultaneously. The signal probability density functions (PDFs) 
are represented by Gaussians for both AE and Mbc. The means and widths are verified 
using B^ — s> D^Ti^ and B J/ipK* decays. Additional bifurcated Gaussians (Gaussians 
with different widths on either side of the mean) are used to model the tails in the AE 
distribution of (pK* channels. The continuum PDF for Mbc {^E) is determined from the 
events outside of AE (Mbc) signal window. The continuum PDFs for Mbc and AE are 
represented by an empirical background function introduced by ARGUS [10| and a linear 
function, respectively. The number of signal and background are floated in the fit while 
other PDF parameters are fixed. The measured branching fractions {B) are summarized in 
Table m The distributions of AE and Mbc for the four measured modes are shown in Fig. Q 

The systematic errors in the signal yields are estimated by varying each fixed PDF pa- 
rameter by ±lcr of its nominal value. Conservatively, the change in the signal yield from 
each variation is added in quadrature. The systematic errors in the efficiency are due to 
uncertainties in track finding (1% per track), particle identification (2%), and vr" finding 
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FIG. 1: Distributions of AE (Mbc) with fit results for the events in the Mj,c (AE) signal window. 
The continuum background component is shown by dashed curves. 



TABLE I: Signal yields (Ng) obtained by fits after background subtraction, total efficiency (e), 
statistical significance (S = ^y2ln[C{Ns)/C{0)]), and measured branching fraction (B). The inter- 
mediate branching fractions are taken from 



Mode 








S 


i3(10-6) 


B+ ^ 






16.9 


16.5 9.4 ± 1.1 ±0.7 


B^^ 




35.6+«-^ 


4.6 


8.7 


9.0+f^ ±0.7 


B^^ 


(l)K*^ 


58.5t^:l 


6.9 


11.3 


^^■^-1.5 -0.8 


B+ ^ 
K*^ 


(j)K*+ 
^ ^ K+vrO 




1.4 
2.1 


4.9 
2.8 
4.0 


7+2.1 +0.Y 
"•'-1.9 -1.0 
f. Q+3.8 +0.9 
'J->^-3.2 -1.0 

c;+2.6 +0.6 
"•"3-2.3 -0.9 



(4%), and to the uncertainty in B{(j) K^K ) (1.4%). The estimated contaminations of 
B foK^*^ and B (pKn are included as an uncertainty in the background. For the 
B — >■ (pK* modes, an additional systematic error in the efficiency due to the uncertainty in 
the polarization together with the uncertainty in the slow pion detection efficiency (l%-4%) 
is included. 

For the self-tagging modes B~^ 4>K~^, B^ (f)K*^{K~^ 7r~), and B~^ (j)K*~^ we have 
studied the direct CP asymmetries Acp = ; where B (B) is B^ or B^ (5° 

or B~) and / is one of the self-tagged (pK^*^ final states. The values of Acp for B^ — >■ 
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(t)K+, 5° (j)K*°{K+7r-), and B+ (f)K*+ are: 0.01 ± 0.12 ± 0.05, 0.07 ± 0.15l|5:Si, and 
—0.13 ± 0.29lo;?^, respectively. These correspond to 90% confidence level limits of —0.20 < 
Acpi(l>K+) < 0.22, -0.18 < Acp{(pK*%K+7t-)) < 0.33, and -0.64 < Acp{^K*+) < 0.36, 
respectively. The systematic error includes the uncertainties in signal extraction (2%) and 
detector induced bias (1 — 6%), which has been studied using large samples of inclusive 
charged kaon and pion tracks, high momentum —>■ K~tt^ and — ^ ir~7r+7r+ decays, 
and 5- meson decays to the channels J /ipK'^*\ D 7r+ and D*~t[^ . The systematic errors due 
to background from B foK^*^ and non-i^'* background in i? ^ (pK* channels are also 
included. 

The decay angles of a 5-meson, to the two vector mesons and K*^, as defined in the 
transversity basis , are shown in Fig. |21 The x-y plane is defined by the K*'^ daughters 
and the x axis is in the direction of the 0-meson. The y axis is perpendicular to the x axis 
and is on the same side as the kaon from decay. The z axis is perpendicular to the x-y 
plane according to the right-hand rule. Otr (0tr) is the polar (azimuthal) angle with respect 
to the z-axis of the from decay in the rest frame. 6k* is defined above. 

The distribution of decays in the three angles 1^1, 9k*, Otr, and 0tr is 



tr, COs6'tr, COs6'i^.) 9 



COS COS 9k* 



327r 



|A_Lp2cos^ 9tr sin^ 9k* 



)tr sin^ 9k* 



+ |y4||p2sin^6'tr sin^ 

+ |v4op4 Sin^ 6'tr COS^ 0tr cos^ 9k* 

+ y2Re(A|Ao) sin^ 9tr sin 20tr sin 29k* 

— riV21m{AQA±) sin 29tr cos 0tr sin 29k* 

— 2r]lm(At,A±) sin 29tr sin 0tr sin^ 9k*] , 



(1) 



where Aq, Ay, and A±_ are the complex amplitudes of the three helicity states in the transver- 
sity basis with the normalization condition |y4oP + |A||p + |A_|_p = 1, and t] = +1 (—1) for S° 
(5°). Aq denotes the longitudinal polarization of the — > K^K~ system and A±_ {A\\) is the 
transverse polarization along the 2-axis (y-axis). The value of \ A\ P (|y4oP+|A||p = 1 — 
is the CP-odd (CP-even) fraction in the decay B (pK* jl3|. The imaginary phases of 
the amplitudes are sensitive to final state interactions (FSI). The presence of FSI results in 
phases that are not either or ±7r. 




FIG. 2: The definition of decay angles in i? — > (pK*^ decay. 



The complex amplitudes are determined by performing an unbinned maximum likelihood 
fit Q with po (l)K*\K+Ti-) candidates in the signal window. The combined likelihood 
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is given by 



N 

L = '[[e{eK*,6tr,4>tr)[UK*0 " r(6'i^. , 6'^, 0tr) 
i 

+ fqg ■ RqqidK* , Otv, 0tr) 

+ fxKK'O ■ RkKK* o(^i^.,^tr,0tr)] , (2) 

where V is the angular distribution function (ADF) given by Eq. HI and Rqq and Rkkk*° 
are the ADFs for continuum and B — > KKK*^ background, respectively. The value of 77 is 
determined from the charge of the kaon in K*^ decay, Rqq is determined from sideband data 
and Rkkk*o is assumed to be flat. The detection efficiency function (e) is determined by 
MC. The fractions of (f)K*^ ifcpK'o), ifgq) and B KKK*^ {fKKK*o) are parameterized 
as a function of AE and M^c- The value of arg(y4o) is set to zero and l^yP is calculated 
from the normalization constraint in the fit. Four parameters (l^oP? arg(A||), and 

arg(y4^)) are left free to be determined from the fit. 

Figure 121 shows projections for each of the three angles together with results from the fit. 
The amplitudes obtained from the fit are |AoP = 0.43±0.09±0.04, \A±\^ = 0.41±0.10±0.04, 
arg(A||) = -2.57 ± 0.39 ± 0.09, and aig{A±) = 0.48 ± 0.32 ± 0.06, where the first errors are 
statistical and the second errors are systematic. The systematic uncertainties include the 
slow pion detection efficiency (3 — 6%), the background from higher K* states (6 — 9%), 
and the B foK* background (1%). The systematic uncertainty due to the angular 
resolution is estimated by MC simulation and found to be less than 1%. Uncertainties due 
to the background PDFs, the signal yields, and the modeling of efficiency function (e) are 
estimated to be 1 — 3%. 




FIG. 3: Projections of the transversity angles with results of the fit superimposed. The points with 
error bars show the efficiency corrected data after subtraction of continuum and B K^K~K* 
background. The x^/n.d.f. for the projection of cos 9k*, cos^tr and (ptr are given by 1.05, 0.90 and 
0.46, respectively. 

In summary, we measure the branching fractions of four B — > (j)K^*'^ decay modes. The 
value of B{B^ — > (j)K~^) is in good agreement with, and supersedes, previously reported Belle 
measurements 0, ■ Our branching fraction results are in agreement with measurements 
by BaBar and CLEO tlSi], and the predictions by PQCD 5]. The measured direct CP 
asymmetries in these modes are consistent with zero. The decay amplitudes for B^ (f)K*^ 
are determined by fitting the angular distributions in the transversity basis. The longitudinal 
polarization fraction {/^{(pK*'^)) reported by BaBar agrees with our measured value 
of |v4op. The measured value of \A±\'^ shows that both CP-odd (|v4_Lp) and CP-even 
(l^oP+l^ll P) components are present in ^i^* decays, in constrast to the case of P — > J/ipK*, 
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which is dominantly CP-even Our data also yield a good fit when the phases of A± and 
A\\ are constrained to zero and — vr, indicating that our data cannot distinguish the presence 
of final state interactions. 
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